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ABSTRACT 



In direct sequence-code division multiple accesses systems (DS-CDMA), the signals are transmitted over 
multipath channels that introduce fading and shadowing. Multipath fading and shadowing along with 
multiple access interference and inter-symbol interference degrades the system performance. Further, 
experimental results have confirmed the presence of impulsive noise in wireless mobile communication 
channels. Hence, this paper presents a technique for multiuser detection in DS-CDMA systems over 
generalized-K (GK) fading channels in presence of impulsive noise. An approximate closed -form expression 
for average probability of error of an M-decorrelator is derived, for the demodulation of binary phase shift 
keying (BPSK) signals transmitted over the channels with simultaneous presence of fading and shadowing. 
Maximal ratio combining (MRC) receive diversity is also incorporated to mitigate the effects of fading and 
shadowing. Performance of a new M-estimator based detector is also studied by evaluating error rate with the 
derived expression. Simulation results reveal that the proposed M-estimator based detector performs better in 
the presence of fading, shadowing and heavy-tailed impulsive noise. 

Keywords: Diversity combining, fading channel;GK distribution; impulsive noise; multiuser detection; 

probability of error, shadowing. 



I. INTRODUCTION 

Recent research has explored the potential benefits of 
multiuser detection for code division multiple access 
(CDMA) communication systems with present 
multiple access interference (MAI) [1]. These optimal 
multiuser detectors have led to the developments of 
the various linear multiuser detectors with Gaussian 
noise though various experimental results have 
confirmed that many realistic channels are impulsive 
in nature [2], [3]. Lately, the problem of robust 
multiuser detection in non-Gaussian channels has 
been addressed in the literature [4], [5], [6], which 
were developed based on the Huber, Hampel, and a 
new M-estimator (modified Hampel), respectively. 
Since CDMA transmissions are frequently made over 
channels that exhibit fading, it is of interest to design 
receivers that take this fading behavior of the channel 
in to account [3]. Robust multiuser detection for 
synchronous DS-CDMA system with maximal ratio 



combiner (MRC) receive diversity over Nakagami-m 
fading channels is presented in [7] by assuming that 
the modulation is binary PSK (BPSK). But, the 
simultaneous presence of multipath fading and 
shadowing leads to worsening the wireless channels 
[8]. Recently, the performance of M-decorrelator in 
simultaneous presence of fading and shadowing 
with impulsive noise is presented in [9]. Hence, in 
this paper, the problem of multiuser detection for 
DS-CDMA system over generalized-K (GK) fading 
channels with impulsive noise is considered. D- 
branch MRC receive diversity is incorporated to 
mitigate the effects of multipath fading. 

DS-CDMA system over multipath fading channel in 
non-Gaussian impulsive noise environment is 
presented in Section II. Section III presents the GK 
fading channel with MRC receive diversity. 
Influence function of the proposed M-estimator is 
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presented in Section IV and an approximate closed- 
form expression for average probability of error of an 
M-decorrelator is also derived. Section V presents 
some simulation results. Section VI concludes the 
paper. 

II. SYSTEM MODEL 

An L-user CDMA system, where each user transmits 
information by modulating a signature sequence 
over a single-path Nakagami-m fading channel, is 
considered in this paper. The received signal over 
one symbol duration can be modeled as [3] 



L M-l 

r(t) = 9U i Z Z (Oe 
1=1 i=0 



A(t) 



s^t-iT.-TM + nit) (1) 



where denotes the real part, M is the number of 

data symbols per user in the data frame of interest, T s 
is the symbol interval, <Z/(f) is the time-varying 

fading gain of the 1 th user's channel, 0/0) is the time- 
varying phase of the 1 th user's channel, (i) is the i th 

bit of the 1 th user, Sj (t ) is the normalized signaling 
waveform of the 1 th user and n(f) is assumed as a 

zero-mean complex non-Gaussian noise. It is 
assumed that the signaling constellation consists of 
non-orthogonal signals given by [3] 

J2 ( ) 7(^+3) for ts[0,T s ] 
s t (t) = WT /v (2) 

0 for m[0,T s ] 

where j = -J— 1 , 0 l is the phase of the 1 th user relative 
to some reference, co c '\s the common carrier 
frequency, and the spreading waveforms, a^t), are 
of the form 

N i 



d\(t) = X a n u T 

n=\ 



where {a l n } n=1 is a signature sequence of +ls and -Is 
assigned to the 1 th user user, and u T (t) is a unit- 

c 

amplitude pulse of duration T c with T s = NT C . At the 

receiver, the received signal r(t) is first filtered by a 
chip-matched filter and then sampled at the chip rate, 
1 / T c . The resulting discrete-time signal sample 

corresponding to the n th chip of the i th symbol is given 
by [3] 

t iTMn+DT. 

rJ<»Jt 



2 iT s +(n+l)T c 

I r(t)e 



T c iT s +nT c 



dt,n = l...N. (4) 



For synchronous case (i.e., tj_ = r 2 = ... = r/ = 0 ), 

assuming that the fading process for each user varies 
at a slower rate that the magnitude and phase can 
taken to be constant over the duration of a bit, (4) can 
be expressed in matrix notation as [3] 

r (i) = A0(i) + w(i) (5) 

where 

T 



r(i) □ \j\(i\...,r N (i) 
w(i) □ \_wi(i),...,w N (i) 

0{i) □ ^[h(i)gl(0^b L (0g L (i)] 



T 



(6) 
(7) 
(8) 



with w n (i) is a sequence of independent and 

identically distributed (i.i.d.) complex random 
variables whose in-phase and quadrature 
components are independent non-Gaussian random 
variables with probability density function (PDF) of 
this noise model has the form 



f = (l- f )X(0, v 2 ) + f X(0, kv 2 ) 



(9) 



with u>0, 0<s<\, and k > 1 . Here K(0,u 2 ) 
represents the nominal background noise and the 
X(0,/tl> 2 ) represents an impulsive component, with s 

representing the probability that impulses occur. It is 
assumed that the 1 th user employs binary phase shift 
keying (BPSK) modulation to transmit the data bits 
bi g [-1,1] with equal probability and a symbol rate 

1/Ts. 

III. GK FADING CHANNEL WITH MRC 
DIVERSITY 

It is assumed that the signal of each user arrives at 
the receiver through an independent, single-path 
fading channel. For the shadowed fading channels, 
a l (i) are i.i.d. random variables with GK distribution 

given by [8] 



T(m)T(ju) 



v 



\mju 



m+/j m+ju 



K 



o J 



m-/j 



V 



Imju 



n 



o 



J 



(10) 

where, m is the Nakagami fading parameter that 
determines the severity of the fading, ju represents 

the shadowing levels, Q 0 is the average SNR in a 

shadowed fading channel, Kz(-) is the modified 

Bessel function and T() is the Gamma function [10]. 

Assuming that the fading is mitigated through D- 
branch MRC, the output of maximal ratio combiner 
can be written as [8, 11] 
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D 

R = 

7=1 



(11) 



where a are i.i.d. GK distributed random variables 
each having a PDF of the form (10). The PDF of R, 
with multipath fading and shadowing from branch to 
branch are distinct, is given by [8] 



Pr(t) = 



T(m m )T(ji m ) 
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where 



m m = Dm + (D - 1) 



and 
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Fig. 1. Influence functions of (a) the linear 
decorrelating detector, (b) Huber estimator, (c) 
Hampel estimator, and (d) the proposed estimator. 

IV. MULTI TAPER SPECTRUM ESTIMATION 
METHOD 

In this section, the proposed M-estimator is 
presented and the average probability of error is 
derived for a single-path shadowed fading channel. 
A) M-estimator 

An M-estimator is, a generalization of usual 
maximum likelihood estimates, used to estimate the 

unknown parameters ^,#2 (where 6=Ab ) by 

minimizing a sum of function /?(•) of the residuals 
[4] 

N I L 

9 = arg min 




Z s\e x 
i=i J 



where p is a symmetric, positive-definite function 

with a unique minimum at zero, and is chosen to be 
less increasing than square and N is the processing 
gain. An influence fun- ction y/{-) = p '(•) is proposed 

(as shown in Fig. 1 along with other influence 
functions), such that it yields a solution that is not 
sensitive to outlying measurements, as [6] 



Vpro (*) = 



x 

asgn(x) 
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v 



X 



2^ 



b* 



for \x\<a 

for a<x<b (16) 



for x >b 



J 



where the choice of the constants a and b depends on 

the robustness measures. 

B) Average probability of error 

The asymptotic probability of error for the class of 
decorrelating detectors, for large processing gain N, 
is given by [4] 

Pj-Pr(4<0|^>0) 

f \ 



= Q 



V 



A l 




R 1 








// J 



(17) 



where Q(x) is the Gaussian Q-f unction [11], 

r 2 

2 )¥ (u)f(u)du 

y = 2 

[\y/\u)f(u)du\ 



(18) 



and R = S S with S is an NxL matrix of columns ai. 
Average probability of error can be obtained by 
averaging the conditional probability of error (17) 
over the PDF, (12), of maximal ratio combiner output 
as 
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Using the 



well known upper-bound approximation, called 
Chernoff bound, to Q(x), given by 



1 -x 2 /2 
Q(x)<-e 



(20) 



(15) the integral of (19) can be expressed as 
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Now, by using [Eq. 6.631.3, 10] to evaluate the 
integral (21), the average probability of error can be 
derived as [9] 
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Fig. 2. Average probability of error versus SNR for 
user 1 for linear multiuser detector (Least Squares), 
minimax detector with Huber, Hampel and 
proposed M-estimator in synchronous CDMA 
channel with impulse noise, N = 31, s = 0. 



V. SIMULATION RESULTS 

In this section, the performance of M-decorrelator is 
presented by computing (22) for different values of 
order of diversity. It is assumed that m = 1.4, |a = 2 
and CI o= 10 dB in the simulations. 

Performance of M-decorrelator with different 
influence functions is shown in Fig. 2, Fig. 3 and Fig. 
4. In Fig. 2, the average probability of error versus 
the signal-to-noise ratio (SNR) corresponding to the 
user 1 under perfect power control of a synchronous 
DS-CDMA system with six users (L = 6) and a 
processing gain, N = 31 is plotted for Gaussian noise 
( s = 0). In Fig. 3, the average probability of error is 
plotted for moderate impulsive noise (s= 0.01). 
Similarly, in Fig. 4, the average probability of error is 
plotted for highly impulsive noise ( s = 0.1). 
Simulation results show that the proposed M- 
estimator based detector performs well compared to 
linear multiuser detector, minimax detector with 
Huber and Hampel estimator based detectors. 



s = 0.01. k = 100. L = 6. Q = 10 dB 




0 5 10 15 20 25 30 35 

SNR in dB 



Fig. 3. Average probability of error versus SNR for 
user 1 for linear multiuser detector (Least Squares), 
minimax detector with Huber, Hampel and 
proposed M-estimator in synchronous CDMA 
channel with impulse noise, N = 31, e=0.01. 



LIRA | 2014 | Volume 1 | Issue 4 



| 171 



International Journal of Research and Applications (Oct - Dec © 2014 Transactions) 
International Conference on Emerging Trends in Electronics & Telecommunications (ICETET-14) 



6 - O.I. k - 100. L - 6. O 10 dU 




SNR in dH 

Fig. 4. Average probability of error versus SNR for 
user 1 for linear multiuser detector (Least Squares), 
minimax detector with Huber, Hampel and proposed 
M-estimator in asynchronous CDMA channel with 
impulse noise, N = 31, e = 0.1. 

VI. CONCLUDING REMARKS 

Multiuser detection for DS-CDMA systems over GK 
fading channels under impulsive noise environment 
is presented. An approximate closed-form expression 
for average probability of error of the M-decorrelator 
to detect BPSK signals is derived by incorporating 
MRC receive diversity. An M-estimator based 
multiuser detector is proposed and its performance is 
analyzed by computing average probability of error 
using the expression derived. Simulation results 
show that the proposed multiuser detector offers 
significant performance gain over the linear 
multiuser detector and the minimax decorrelating 
detector with Huber and Hampel M-estimators, in 
heavy-tailed impulsive noise. 

REFERENCES 

[1] S. Verdu, Multiuser Detection, Cambridge Univ. 
Press, Cambridge, U.K., 1998. 

[2] A.M. Zoubir, V. Koivunen, Y. Chakhchoukh, and 
M. Muma, "Robust Estimation in Signal 
Processing: A Tutorial-Style Treatment of 
Fundamental Concepts," Signal Processing 
Magazine, IEEE, vol.29, no.4, pp.61-80, Jul. 2012. 

[3] H. V. Poor, and M. Tanda, "Multiuser Detection 
in flat fading non-Gaussian channels/' IEEE 
Trans. Commun., vol.50, no. 11, pp. 1769-1777, Nov. 
2002 



[4] X. Wang and H. V. Poor, "Robust multiuser 
detection in non-Gaussian channels," IEEE Trans. 
Signal Process., vol. 47, no. 2, pp. 289-305, Feb. 
1999. 

[5] T. Anil Kumr, K. Deergha Rao, and M.N.S. 
Swamy, "A robust technique for multiuser 
detection in non-Gaussian channels,"in Proc. 47 th 
Midwest Symposium on Circuits and Systems, vol. 3, 
pp. iii-247-50, vol. 3, 25-28, Jul. 2004. 

[6] T. Anil Kumar, and K. Deergha Rao, "Improved 
robust techniques for multiuser detection in non- 
Gaussian channels", Circuits Systems and Signal 
Processing J., vol. 25, no. 4, 2006. 

[7] V. Srinivasa Rao, P. Vinay Kumar, S. Balaji, 
Habibulla Khan, and T. Anil Kumar, "Robust 
multiuser detection in synchronous DS-CDMA 
system with MRC receive diversity over 
Nakagami-m fading channel", Advanced 
Engineering Forum Smart Technologies for 
Communications, vol. 4, pp. 43-50, Jun. 2012. 

[8] P. M. Shankar, "Maximal Ratio Combining (MRC) 
in Shadowed Fading Channels in Presence of 
Shadowed Fading Cochannel Interference(CCI)," 
Wireless Pers Commun J., vol. 68, no. 1, pp. 15-25, 
2013. 

[10] V. Srinivasa Rao, Habibulla Khan and T. Anil 
Kumar, "Multiuser detection in shadowed 
fading channels with impulsive noise", Proc. of 
Int. Conf. on Emerging Trends in Engg. and Tech 
(ICETET 2014)., vol 1. , pp. 36-40, 29-30 Sep. 2014. 

[11] I. S. Gradshteyn, I. M. Ryzhik, Table of Integrals, 
Series, and Products, Academic Press, Boston, 
MA, 2007. 

[12] P. Sharma, "Selection of diversity and 
modulation parameters for Nakagami fading 
channels to jointly satisfy outage and bit error 
requirements," IEEE Transactions on Wireless 
Communications, vol.5, no. 6, pp. 1279-1283, June 
2006. 

[13] J.G. Proakis, Digital Communications, Fourth 
Edition, McGraw Hill, 2001. 



LIRA | 2014 | Volume 1 | Issue 4 | 172 



